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EXECUTIVE  SUMMARY 


A  number  of  papers  have  appeared  in  the  open  literature  describing  either  experimental  or 
theoretical  work  using  laser  induced  filaments  to  guide  energy  in  the  form  of  microwave  radiation.  The 
idea  is  to  minimize  the  natural  divergence  of  a  radiation  source  by  confining  the  energy  flow  in  a 
specified  direction  with  the  aide  of  conducting  plasma  filaments  which  would  maintain  a  high  energy 
density  over  a  large  distance.  The  filaments  act  either  as  an  open  waveguide  (Goubau  line)  or  they  can  be 
arranged  to  form  a  hollow  plasma  tube.  The  theoretical  work  published  to  date  in  the  open  literature 
assumes  steady-state  plasma  with  constant  conductivity.  In  reality,  the  plasma  decays  within  tens  of 
nanoseconds  due  to  recombination  and  electron  attachment.  Thus,  the  correct  theoretical  treatment  of  the 
problem  should  take  into  account  the  transient  effects  involved  in  both  the  launching  of  the 
electromagnetic  wave  and  the  simultaneous  loss  of  plasma  conductivity.  The  transient  launching  problem 
is  considered  here. 

An  analytic  model  to  describe  the  transient  propagation  of  an  electromagnetic  wave  inside  a 
hollow  constant  conductivity  plasma  annular  tube  is  formulated.  The  model  neglects  external  radiative 
emission  since  the  plasma  tube  wall  is  assumed  to  be  several  skin  depths  in  thickness.  Approximate 
analytic  time-dependent  solutions  for  mode  propagation  are  obtained  by  a  combination  of  integral 
transform  techniques  and  the  method  of  eigenfunction  expansion  with  homogeneous  boundary  conditions. 
The  model  provides  insight  into  wave  attenuation,  the  dispersion  relations,  waveguide  excitation  and 
mode  coupling  due  to  finite  wall  plasma  conductivity.  The  model  is  developed  by  first  considering  the 
transient  propagation  characteristics  of  ideal  cylindrical  waveguides  with  perfectly  conducting  walls. 

This  preliminary  analysis  provides  a  consistent  notation  for  the  closed  plasma  waveguide  and  suggests 
possible  launching  modes  with  minimal  loss.  The  initial  condition  imposed  on  the  annular  plasma 
waveguide  is  fixed  by  exciting  a  specified  mode  in  the  ideal  cylindrical  waveguide  which  is  then 
transferred  into  the  plasma  tube. 

The  propagation  problem  with  Joule  heating  within  the  plasma  tube  wall  is  analyzed  numerically 
using  a  finite-element  technique  and  some  preliminary  results  are  given. 
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I.  INTRODUCTION 

A  number  of  experimental  and  theoretical  papers  have  appeared  in  the  open  literature  describing 
the  use  of  laser- induced  filaments  to  guide  energy  in  the  form  of  microwave  radiation1'5.  The  natural 
divergence  of  a  radiation  source  can  be  reduced  if  the  energy  flow  is  confined  within  a  limited  spatial 
region  by  walls  (closed  waveguide)  or  a  guiding  structure  such  as  a  parallel  transmission  line.  The 
simplest  open-guide  structure  is  a  single  dielectric  coated  wire  of finite  conductivity  known  as  a  Goubau 
line6.  The  electromagnetic  energy  is  transported  along  this  open  line  as  an  attached  surface  wave.  The 
propagation  of  the  surface  wave  depends  upon  the  ability  of  the  electromagnetic  field  to  penetrate  the 
conductor’s  interior.  A  perfect  single  conductor  will  not  support  such  a  guided  wave.  The  propagation 
characteristics  of  the  guided  wave  are  determined  by  the  wire  conductivity,  the  wire  radius,  and  the 
radiation  frequency.  The  guiding  performance  is  enhanced  by  coating  the  wire  with  a  dielectric.  The 
calculation  of  the  propagation  characteristics  is  based  upon  tangential  continuity  of  the  electric  and 
magnetic  field  components  across  the  dielectric -wire  surface  boundary  and  utilizes  a  steady-state 
sinusoidal  approximation.  The  standard  analysis  neglects  how  the  wave  was  launched. 

The  laser-induced  plasma  filaments  behave  as  ohmic  resistors  over  distances  on  the  order  of  tens 
of  centimeters7.  By  treating  a  single  filament  as  a  time-dependent  resistance,  we  have  determined  the 
electron  density,  the  kinetic  parameters  for  electron  attachment  and  recombination,  and  an  effective 
conductivity  of  350  S/m.  In  order  to  arrive  at  this  conductivity,  we  assumed  that  the  filament  was  a 
cylinder  with  a  diameter  on  the  order  of  100  pm  and  that  the  plasma  electrons  had  a  collision  frequency  of 
1012  s'1.  Later  experiments  conducted  at  NRL  obtained  the  conductivity  by  a  technique  which  does  not 
require  knowledge  of  the  electron  collision  frequency.  The  initial  filament  conductivity  is  50  S/m  and 
decays  to  4-5  S/m  over  a  few  hundred  nanoseconds. 

The  laser  is  a  conventional  chirped  pulse  amplified  Ti:  Sapphire  laser.  A  KML  oscillator  is 
stretched  and  amplified  to  about  1.5  mJ  in  a  Coherent  Spitfire  regenerative  amplifier  at  1  kHz.  These 
pulses  are  then  amplified  by  a  Coherent  custom  built  frequency  doubled  Nd:YAG  laser  in  a  four  pass 
amplifier  to  about  1.2  Joules  in  ~260  picosec  at  10  Hz  repetition  rate.  The  stretched  pulses  are  expanded 
to  40  mm  diameter  and  compressed  with  dual  gold  coated  gratings  resulting  in  40  femtosec,  850  mJ 
pulses.  A  simple  two  mirror  telescope  is  used  to  control  the  spot  size  down  stream  from  the  laser  and  the 
position  of  filament  formation  is  controlled  by  adjusting  the  grating  separation.  Typically,  spot  sizes  on 
the  order  of  1  -  1.5  cm  diameter  are  used  in  the  experiments  and  contain  70-80  filaments.  Although  the 
filament  distribution  can  change  from  day  to  day,  normally  once  the  system  has  stabilized  (after  40 
minutes)  the  distribution  remains  the  same  for  the  rest  of  the  day. 

In  any  experiment  attempting  to  use  the  laser-induced  filaments  as  a  guiding  structure,  the 
transient  behavior  of  the  plasma  and  the  electromagnetic  interaction  within  the  particular  experimental 
configuration  must  be  understood  in  order  to  correctly  interpret  the  results.  The  NRL  femtosecond  laser 
is  capable  of  consistently  producing  75-80  filaments  within  a  one  centimeter  diameter  disk  at  ten  hertz 
frequency  over  several  meters.  The  filaments  are  “produced”  on  a  time  scale  of  a  few  nanoseconds  over 
the  experimental  geometry.  The  electromagnetic  field  measurements  are  on  the  same  time  scale.  Thus, 
all  of  our  observations  involve  transient  processes  and  the  standard  steady-state  approximation  of 
sinusoidal  waves  is  inapplicable  in  theoretical  or  experimental  work  with  short-lived  laser-induced 
filaments.  This  report  documents  our  current  progress  in  developing  such  transient  models. 

An  alternative  experimental  approach  which  can  avoid  the  transient  plasma  lifetime  to  some 
extent  is  to  create  the  plasma  in  air  by  focusing  high-energy  laser  beams.  There  are  several  published 
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papers  that  use  this  technique8"12.  Plasmas  created  in  this  manner  do  not  contain  the  laser-generated 
filaments  produced  by  self-focusing  and  require  enormous  energy.  The  electron  density  is  lower  and 
creation  of  the  plasma  may  require  seeding  the  atmosphere  with  a  hydrocarbon.  Microwave  radiation  has 
been  guided  over  a  distance  of  60  m  at  35.3  GHz  frequency  by  a  closed  waveguide  created  in  this 
fashion8'9.  However,  we  believe  that  the  energy  cost  (50  J)  and  the  complexity  of  the  KrF  excimer  laser 
used  in  the  experiment  as  well  as  the  required  optics  would  be  offset  by  simply  using  a  larger  microwave 
source. 


In  order  to  develop  a  transient  model  of  the  electromagnetic  interaction  with  the  plasma,  an 
experimental  configuration  must  be  selected  for  the  analysis.  Since  we  are  interested  in  guiding  the 
energy,  a  simple  circular  metal  waveguide  coupled  to  a  hollow  plasma  tube  whose  walls  are  several  skin 
depths  in  thickness  was  selected  as  a  preliminary  base  line  model.  The  metal  waveguide  provides  an 
initial  condition  for  the  propagation  of  the  wave  into  the  plasma  tube.  If  the  experimental  geometry  were 
changed  to  a  circular  horn  for  example,  the  same  mathematical  technique  described  here  should  still  be 
applicable.  The  basis  sets  used  in  the  model  would  be  modified  to  the  appropriate  geometry  and 
coordinate  system. 

Section  II  develops  the  notation,  dispersion  relations,  and  description  of  the  modes  and  cut  off 
frequencies  for  the  ideal  circular  waveguide.  This  notation  is  utilized  in  Section  III  for  the  development 
of  the  time-dependent  theory.  If  the  reader  is  familiar  with  these  concepts  he  should  be  able  to  skip 
Section  II  and  proceed  to  the  transient  analysis  in  Section  III. 


II.  STEADY-STATE  WAVEGUIDE  THEORY 

The  propagation  of  electromagnetic  radiation  in  a  lossy  waveguide  is  governed  by  Maxwell’s 
equations.  The  analytic  solution  of  the  Maxwell  equations  provides  the  axial  propagation  constant  for  the 
radiation  field,  the  eigenvalues  and  eigenfunctions  associated  with  a  particular  mode  for  the  specified 
geometry  (waveguide  cross  section),  and  the  waveguide  dispersion  relation  as  a  function  of  geometry  and 
frequency.  The  attenuation  due  to  finite-wall  conductivity  is  often  treated  as  a  perturbation  of  the  perfect 
wall  conductor  solutions.  Transient  effects  are  customarily  neglected,  since  steady-state  sinusoidal 
solutions  are  most  often  of  practical  engineering  interest.  Losses  due  to  the  dielectric  (air)  within  the 
guide  are  usually  negligible. 

The  standard  engineering  procedure  in  solving  the  waveguide  problem  for  a  sinusoidal  input 
source  is  to  transform  the  Maxwell  equations  to  a  coordinate  system  appropriate  to  the  waveguide 
geometry  and  then  assume  that  the  waveguide  walls  are  perfect  conductors14.  The  assumption  of  perfect 
wall  conductivity  simplifies  the  imposition  of  the  boundary  conditions.  The  tangential  component  of  the 
electric  field  and  the  normal  component  of  the  magnetic  field  intensity  must  vanish  at  the  waveguide  wall. 
These  boundary  conditions  impose  the  mode  structure  and  “restrict”  the  propagation  of  the  radiation  field 
to  the  interior  of  the  guide.  In  reality  the  electromagnetic  wave  penetrates  the  wall,  but  is  exponentially 

attenuated  over  a  radial  distance  of  a  few  skin  depths.  The  skin  depth  S(f,  a)  is  given  by 


S(f,a)  = 


1 


p07rfa 


(1.1) 
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f  denotes  the  radiation  frequency  and  cr  the  wall  conductivity.  Equation  (1.1)  shows  that  at  higher 
frequencies  and/or  wall  conductivities  the  skin  depth  decreases  and  there  is  less  wall  penetration  of  the 
radiation  field.  At  gigahertz  frequencies  the  skin  depth  in  a  good  conductor  such  as  copper  with 
conductivity  «  107  S/m  is  on  the  order  of  microns.  Note  that  the  waveguide  wall  thickness  should  be  at 
least  2-3  times  the  skin  depth  in  order  to  confine  the  radiation  field  and  provide  a  guiding  structure.  If  the 
wall  conductivity  is  too  low,  the  radiation  field  can  “leak”  out  of  the  waveguide  into  the  surrounding  free 
space. 

In  the  analysis  of  guided  waves  along  uniform  systems,  it  is  customary  to  classify  the  wave 
solutions  into  the  following  three  types: 

1 .  Transverse  electromagnetic  waves  (TEM)  These  waves  have  their  electric  and  magnetic  field 
components  entirely  in  the  transverse  plane.  There  is  neither  an  electric  nor  magnetic  field 
component  in  the  propagation  direction.  TEM  waves  propagate  on  multi-conductor  guides  such 
as  transmission  lines  and  are  sometimes  called  principal  waves.  The  TEM  mode  does  not  exist 
inside  a  single,  hollow,  cylindrical  conductor  of  infinite  conductivity  and  will  not  be  considered 
any  further. 

2.  Transverse  magnetic  waves  (TM)  In  a  TM  wave  there  is  an  electric  field  component  (E_ )  along 
the  direction  of  wave  propagation.  The  magnetic  field  component  lies  entirely  within  the 
transverse  plane.  TM  waves  may  also  have  transverse  electric  components,  e.g.,  Ex,E  or  E, . 

These  waves  are  sometimes  called  E  waves  in  the  literature. 

3.  Transverse  electric  waves  (TE)  In  a  TE  wave  there  is  an  axial  component  of  the  magnetic  field 
(Hz)  along  the  propagation  direction.  The  transverse  plane  may  have  both  electric  and  magnetic 
field  components.  These  waves  are  also  referred  to  as  H  waves. 

We  will  first  review  some  elementary  results  for  a  hollow,  cylindrical  conductor  of  infinite 
conductivity.  These  results  provide  a  zero-order  solution  for  the  transient  model,  establish  the  notation 
and  give  insight  into  the  mathematical  issues  that  arise  in  developing  the  transient  model  of  the  plasma 
waveguide.  Moreover,  a  cylindrical  metal  waveguide  provides  a  possible  launching  device  for  a  specified 
mode.  The  judicious  selection  of  a  low-loss  launching  mode  may  enhance  the  performance  of  the  plasma 
waveguide. 


Electric  and  magnetic  fields  are  vector  quantities,  which  are  related  in  space  and  time  via 
Maxwell’s  equations 


V  ■  D  =  p 
V-B  =  0 


VxE 


8B_ 

dt 


VxH  =  J  + 


8D 

dt 


D  =  s0E 
B  =  p0H 


(2.1) 


The  standard  mathematical  technique  in  solving  the  steady-state  sinusoidal  waveguide  problem  is  to 
assume  a  priori  that  the  vector  solution  can  be  written  in  the  following  form 


E(x,y,z,t )  =  E(x,  y)exp[j cot  - k z] 
B(x,  y,  z,  t )  =  B(x,  v)exp  [jcot  -  kz] 


(2.2) 
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Traveling  or  standing  waves  in  the  axial  direction  (z)  can  be  described  by  taking  appropriate  linear 
combinations  of  the  exponential  functions.  The  parameter  k  is  the  propagation  constant  or  wave  number 
and  may  be  real  or  complex.  Substitution  of  equation  (2.2)  into  (2.1)  yields  a  two-dimensional  vector 
Helmholtz  equation.  The  dispersion  relation  k(co)  between  the  wave  number  k  and  the  sinusoidal 
frequency  CO  is  obtained  by  solving  this  Helmholtz  equation  with  appropriate  boundary  conditions  for 
TM  or  TE  modes.  In  a  waveguide  with  walls  of  infinite  conductivity  the  boundary  conditions  are 
imposed  at  the  inner  radius  of  the  cylindrical  wall  and  result  in  a  well-posed  eigenvalue  problem  whose 
solution  gives  the  dispersion  relation.  If  the  skin  depth  is  small  in  comparison  with  the  wall  thickness  and 
waveguide  radius,  finite  conductivity  effects  can  be  incorporated  into  the  analysis  by  perturbation  theory 
or  a  variational  calculation.  The  incoiporation  of  finite  wall  conductivity  yields  a  complex  wave  number 


k{(0) 


k(co)  =  j[k0{co)+P{co)\  +  a{co) 


(2.3) 


where  k°(co)  is  the  wave  number  for  perfectly  conducting  walls  and  p(co)  is  the  change  in  the 
propagation  constant  due  to  finite  conductivity.  In  a  good  conductor,  the  /?(&>)  correction  is  usually  small 
and  can  be  neglected.  However,  the  fHco)  term  may  be  significant  near  the  cut-off  frequency  C0cut , 
where  k°  ( cocut )  =  0 .  The  a  ( a)  term  is  the  attenuation  coefficient  due  to  conduction  currents  within  the 

wall  which  result  in  Joule  heating.  The  dispersion  relations  for  a  single,  hollow,  cylindrical  conductor  of 
inner  radius  b  and  infinite  conductivity  are  given  by15 


km  ( ('j ) 


If 3 

2 

( P 

\ 

— 

— 

l  b 

) 

Ke{®)  =  \  - 


If  a?) 

2 

fa.) 

I{cj 

1  b  ) 

(2.4) 


where  P  and  0„„  denote  the  m-th  root  of  the  Bessel  functions  J  „  ( Pm )  =  0  and  J'  (Q„„, )  =  0 

nm  z~-nm  n  \  nm  ’  n  ' 

respectively.  The  integer  n  denotes  the  order  of  the  Bessel  function.  Values  of  these  roots  are  given  in 
Table  I. 


TABLE  I  -Bessel  Function  Zeros 


MODE 

Pnm 

MODE 

Qnm 

TMoi 

2.405 

TE„ 

1.841 

TM„ 

3.832 

te21 

3.054 

tm21 

5.136 

TEoi 

3.832 

TM02 

5.520 

te12 

5.331 
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The  frequency  at  which  k°  (co)  is  zero  is  the  cut-off  frequency  (Ocut  =  2 K  fcut .  Modes  with  frequencies 

below  the  cut-off  frequency  will  not  propagate  in  the  perfect  waveguide  since  the  wave  number  k"  ( co)  is 

purely  imaginary  and  the  electromagnetic  wave  amplitude  is  exponentially  attenuated  with  axial  distance 
z.  This  attenuation  is  not  dissipative  as  is  that  associated  with  wall  resistance  and  conductance  in  systems 
with  propagating  waves.  Below  the  cut-off  frequency,  the  waveguide  is  analogous  to  a  filter  composed  of 
reactive  elements.  The  electromagnetic  energy  is  not  irreversibly  lost  via  Joule  heating  in  the  wall,  but  is 
reflected  back  to  the  source.  The  waveguide  acts  as  a  pure  reactance  to  a  source  driven  below  the  cutoff 
frequency. 

Table  II  provides  the  cut-off  frequencies  fcut  calculated  from  Equation  (2.4)  for  a  1.0-cm  inner 
diameter  hollow  cylindrical  waveguide  of  infinite  conductivity. 


TABLE  11-Cut-off  Frequencies  fCut(b  =  0.5  cm) 


Mode 

fcut  GHz 

Mode 

fcut  GHz 

TMoi 

22.95 

TEn 

17.57 

TMn 

36.57 

te21 

29.14 

tm21 

49.01 

TEoi 

36.57 

tm02 

52.68 

TE12 

50.87 

The  selection  of  a  1.0  cm  diameter  waveguide  is  dictated  by  experimental  limitations  on  the  laser  beam 
size  and  the  number  of  conducting  filaments  that  can  be  generated  within  the  laser  beam  radius.  The  TEn 
mode  (H  wave)  has  the  lowest  cut-off  frequency  and  is  the  fundamental  mode  of  the  circular  waveguide. 
In  a  waveguide  with  perfectly  conducting  walls  the  eigenfunctions  (axial  modes)  are  orthogonal.  The 
proof  of  this  assertion  relies  upon  the  imposed  boundary  conditions  at  the  perfectly  conducting  wall  and 
the  two-dimensional  form  of  Green’s  second  identity.  However,  due  to  skin-depth  effects  in  actual 
conductors,  in  reality  the  modes  are  not  orthogonal  and  may  couple.  The  dispersion  relations  calculated 
from  equation  (2.4)  for  both  the  TM  and  TE  modes  in  a  cylindrical  waveguide  with  b=0.5  cm  are  given  in 
Figures  1  and  2.  In  the  ideal  waveguide,  the  excitation  of  a  single  mode  at  some  excitation  frequency  f 
would  result  in  that  particular  mode  propagating  with  the  corresponding  wave  number 
kjM  TE  (/)  determined  from  the  mode  dispersion  relation.  If  the  excitation  frequency  is  above  the  cut-off 

frequencies  of  available  lower  modes,  these  available  modes  may  be  excited  as  well.  In  general,  the  net 
effect  of  the  coupling  is  an  increase  in  energy  loss  to  the  surrounding  wall. 

As  an  example,  if  one  were  to  excite  the  TE0i  mode,  which  has  a  cut-off  frequency  of  36.57  GHz, 
with  an  appropriate  source  at  f  =  94  GHz,  the  TE0i  mode  would  propagate  within  the  waveguide  and 
have  a  k°E  =1815  m"1  (Figure  2).  The  dispersion  relation  indicates  that  the  TEn,  TE2i  and  TEn  modes 

could  be  excited  at  the  same  frequency  which  lies  above  all  of  their  cut-off  frequencies  (Table  II). 
Coupling  to  these  modes  via  wall  conductivity  will  weakly  attenuate  and  distort  the  input  signal. 
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Figure  1  Dispersion  relation  for  TM  modes  in  a  cylindrical  waveguide  of  radius  0.5  cm. 


Figure  2  Dispersion  relations  for  TE  modes  in  cylindrical  waveguide  of  radius  0.5  cm. 
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The  reciprocal  slope  of  the  dispersion  relations  provide  information  about  the  group  and  phase 
velocities  of  the  wave  in  the  waveguide.  For  propagating  modes  withf  >  fcut ,  the  phase  and  group 
velocities  are  given  by16 


*  phase 


(f)= 


2nf 


TM.TE 


l- 


(?  \ 
cut 

V  f  J 


group 


(f)  = 


^  3co  ^ 


5k 


o 


/TM.TE 


(f  > 
=  cjl-  -^- 

l  f  J 


(2.5) 


The  electromagnetic  wave  attenuation  for  a  hollow  cylindrical  conductor  of  finite  conductivity 
can  be  calculated  from  the  power  loss  to  the  conductor  wall  by  perturbation  theory  using  the  perfectly- 
conducting  eigenfunctions  as  zero-order  solutions.  The  validity  of  the  calculation  depends  upon  the  skin 
depth  being  small  in  comparison  to  the  wall  thickness.  The  attenuation  coefficient  for  the  TE  and  TM 
modes  in  a  cylindrical  waveguide  is  given  by17 


aTE(n,m,f ) 


aTM(n,m,f) 


0 

1 

{ fa, ) 

bqTo 

V 

f  17  Y 

Y  J  cut 

1/ J 

If) 

yjnfM o 

1 

1 

.  1- 

(  fa,  ) 

i 

l  Z  J 

l^  =  376.73Q 


(2.6) 


These  formulae  are  derived  by  calculating  the  time-averaged  power  flow  into  the  conductor  wall.  The 
power  loss  can  be  related  to  an  effective  surface  current  and  ohmic  dissipation  of  the  electromagnetic 
energy.  In  the  limit  of  f  »  fcut ,  the  attenuation  coefficients  approach  the  following  asymptotic  forms 
for  n  =  0 


a  (Om  /■)=>&  A 

UTE  )  I —  3/2 

ar/yja  J 

„  ,n  ...  ^ 

OL TM  ^5  f  )  / — 


(2.7) 


It  is  to  be  noted  that  the  attenuation  coefficient  for  all  TM  modes  increases  with  the  square  root  of 
frequency.  However,  the  attenuation  coefficient  for  all  TE0m  modes  decreases  as  the  frequency  is 
increased. 
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The  derivation  of  these  attenuation  formulae  relies  upon  the  skin  depth  being  small  in  comparison  to  the 
waveguide  wall  thickness.  The  spatial  variation  of  the  fields  normal  to  the  conducting  surface  is  much 
more  rapid  than  the  field  variations  parallel  to  the  surface.  Figure  3  shows  calculated  skin  depths  as  a 
function  of  frequency  and  conductivity.  In  this  figure  we  have  utilized  the  conductivity  of  copper  (5.8  x 
107  S/m)  as  a  reference  point  and  simply  reduced  the  value  by  a  factor  of  103  or  107.  The  lowest 
conductivity  (5.8  S/m)  roughly  corresponds  to  the  conductivity  of  seawater  and  the  skin  depth  is  on  the 
order  of  millimeters.  Thus  for  materials  of  this  conductivity,  the  waveguide  wall  thickness  should  be 
several  millimeters  thick  to  confine  the  radiation. 


Calculated  attenuation  coefficients  (decibel/m)  are  presented  in  Figures  4  and  5.  These 
calculations  were  performed  for  a  0.5  cm  radius  cylindrical  copper  waveguide  with  a  wall  conductivity 
of  oCu  =5.8x1  ()'  S  /  m  .  Note  the  decrease  in  the  attenuation  of  the  TE0i  mode  (green)  with  increasing 
frequency.  The  attenuation  coefficients  at  a  different  electrical  conductivity  cr  may  be  estimated  from 


either  of  these  figures  by  multiplying  the  y-axis  by  the  quantity 


i 


cr, 


Cu 


However,  the  extrapolation  for 


cr 


attenuation  coefficients  at  very  low  conductivity  is  suspect.  As  the  wall  conductivity  is  reduced,  the  skin 
depth  increases  and  the  assumption  of  a  thin  transitional  layer  with  rapid  decrease  in  the  normal  field 
components  becomes  invalid.  An  alternative  perturbation  method  for  calculating  the  loss  due  to  the 
surrounding  wall  will  be  developed  below. 
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Figure  4  Attenuation  coefficients  (db/m)  for  TM  modes  in  cylindrical  copper  waveguide  of  radius  0.5  cm. 

Note  the  increase  in  attenuation  with  frequency. 


0  20  40  60  80  100  120  140  160 


Frequency  (GHz) 


Figure  5  Attenuation  coefficients  (db/m)  for  TE  modes  in  a  cylindrical  copper  waveguide  of  radius  0.5  cm. 
Note  the  decrease  in  attenuation  with  frequency  for  the  TE0i  mode. 
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In  order  to  develop  a  time-dependent  perturbation  theory  for  the  attenuation  coefficient  in  a 
circular  cylindrical  waveguide,  we  must  first  convert  Maxwell’s  equations  into  time-dependent  wave 
equations.  Taking  the  curl  operator  (  V  x  )  on  the  Maxwell  equations  (2.1)  results  in  the  following  vector 
equations 


V2E  = 


dE 

dt 


1  d2E 

"7- lie 


V2/J  =  /v T^  + 
Ot 


1  82H 


(2.8) 


If  we  consider  TM  or  TE  waves  in  a  cylindrical  waveguide,  the  respective  axial  field  components  obey 
the  following  scalar  equations 


)_d_ 

r  dr 


81E 

dr 


+  - 


1  d2E  d2E 

-  +  - 


dp 


dz 2 


=  °  it)  ■ 


dE  1  d2E 
-  +  ■ 


dt 


dt 2 


]_d_ 
r  dr 


dH. 

dr 


1  d2Hz  d2Hz 
+  — r - - f1 


dp 


dz2 


dH  1  d2H 
=  a//0 — -  +  - 


dt 


dt 2 


(2.9) 


These  equations  are  to  be  solved  for  E _  (r,  p  z,  t)  or  (r,  p  z,  t)  subject  to  the  following  boundary 

conditions  at  the  waveguide  wall  r  =  b 

Ez=  0  ^  =  0  (2.10) 

dr 

These  boundaiy  conditions  correspond  to  a  perfect  conductor  at  r  =  b .  If  the  solutions  are  assumed  to  be 
of  the  form  of  equation  (2.2)  with  perfectly  conducting  walls  [  P(co)  =  0,  a(co)  =  0,  a  =  0  ],  the  axial  field 
components  are  given  by 


Ez(r,(j),z,t)  =  AJn 


H_(r,<j),z,t)  =  AJn 


P  r 

nm 


COS  [lip 

sin(»^) 


exp  (/[ojt-kpz]) 


Q  r 

Z-'tim 


COSpp 

sin(«^) 


exp  (j[cot-kaTEz]) 


(2.11) 


These  solutions  represent  stationary  waves  traveling  in  the  axial  direction.  The  excitation  source  is  at 
minus  infinity.  The  dispersion  relations  (2.4)  are  obtained  directly  from  these  steady-state  solutions. 
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The  transverse  field  components  ( Er ,  ,  1 1  r ,  II  (lt )  may  be  obtained  from  equations  (2.11) 

through  Maxwell’s  equations18 


Er=^ 


5Ez  cofu0  SHz 


_ !_  4. 


Sr  r  Scj) 


E*=i? 


-kxM  gEz 
r  Scj) 


+  CO(U0- 


SH 


Sr 


cos0  SEZ  o  SH 
kte 


r  Scj) 


Sr 


f  \2 

'  (0  ' 


C08n 


SE, 


-  + 


kuTE  SHZ 


Sr  r  Scj) 


In  the  notation  of  reference  18,  the  propagation  constant  y  =  jk°M  TE  . 


(2.12) 


III.  TRANSIENT  ANALYSIS  OF  A  CYLINDRICAL  PLASMA  WAVEGUIDE 

For  simplicity,  we  will  consider  only  n=0  modes  which  implies  that  the  fields  are  independent 
of  (f) .  Since  the  differential  operator  in  (2.9)  is  identical  for  both  axial  field  variables  ( Ez ,  IE  ),  we  may 

represent  the  mathematical  problem  by  a  single  differential  equation  subject  to  appropriate  boundary  and 
initial  conditions 


1  S  f 

r 

V  dr 


r  dr 


dy/\  d  y/ 

j  dz 2 


du z  1  SV 
—  +  ' 


dt 


dt 2 


cr  =  crn 


Y(r,0,t)  =  F(r,t)  ^(r,z,0)  =  0 


r  dy/^ 
dt 


=  0 


Jt= 0 


a  <r  <b  cr  =  0  r  <a 
Y=Y(r,z,t)  0<z<oo 


(3.1) 


Z 


Figure  6  Cylindrical  waveguide  cross-section. 
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Figure  6  shows  a  radial  slice  of  the  cylindrical  waveguide.  The  excitation  source  or  incident  field 
at  z  =  0  is  represented  by  the  function  F(rJ)  =  Fx  (r)S(cot) .  The  plasma  wall  thickness  is  assumed  to 
be  large  enough  so  that  the  tangential  electric  field  is  zero  at  r  =  b  ( z,t)  =  0  ).  The  boundary 
condition  at  r  =  b  for  the  magnetic  field  intensity  (TE  modes)  is  given  by  equation  (2. 10).  The  specified 
initial  conditions  correspond  to  no  field  within  the  waveguide  region  z  >  0  at  /  =  0  . 

In  order  to  solve  equation  (3.1)  we  introduce  the  asymmetric  Fourier  Sine  transform  over  z 
defined  by  Haberman19 


2 

7T 


CO 


z)  sin[A:.z]  dz 


f(r,t,z)  =  |  f(r,  t,k_)  sin  [kzz]  dk . 
0 


(3.2) 


Taking  the  Fourier  Sine  transform  z  — >  k„  of  equation  (3.1)  yields  the  following  inhomogeneous  partial 
differential  equation  for  y/(r,t,kz) 


1  d  f  dy/^ 

r  dr\  dr 


+ — kzFx  (r )S(cot)  =  a(r)ju0  ^  + 
n  ot 


1  d2y 7 


(3.3) 


The  time-dependent  excitation  source  of  frequency  <x>  =  In  f  at  z  =  Onow  appears  explicitly  in  the 
equation.  The  boundary  condition  imposed  upon  the  electric  field  ( i//(b,  t,z)  =  0 )  suggests  that  we  seek 
the  solution  of  equation  (3.3)  for  TM  modes  in  the  form  of  an  eigenfunction  expansion  using  zero-order 
Bessel  functions  of  the  first  kind  J0  (x)  and  time-dependent  coefficients  An  (/,  A:. ) 

oo  oo 

=  F(r,0  =  S(®<)Er»-7«(¥>  W„)  =  0  (3.4) 

n= 1  m= 1 

where  P0n  denotes  the  nth  zero  (Table  I)  of  J0(x) .  The  selected  eigenfunctions  are  finite  at  r  =  0  and 
satisfy  the  imposed  boundary  condition  on  the  electric  field  E,  at  r  =  b  .  The  eigenfunctions  for  TE 
modes  (See  Eqn.  2. 1 1)  would  be  of  the  same  form,  except  that  the  parameter  P0n  would  be  replaced 
by  Q[ht  from  Table  I.  The  excitation  source  F(r,t)  is  also  expanded  in  equation  (3.4)  as  a  Fourier-Bessel 
series  with  an  assumed  time  dependence  of5'(<yf)  .  The  Fourier  amplitudes  Tm  of  the  excitation  source 
are  calculated  from  the  following  integral 

r"  =  4-  j  Nm  =  i  b2j;-(P,J  (3.5) 

Nm  o  b  2 

where  N m  is  a  normalization  factor.  In  general  the  Fourier  amplitudes  may  be  functions  of  time.  We 

shall  assume  that  the  time-dependence  is  factorable  from  the  spatial  dependence  in  the  excitation  source. 
Figure  7  illustrates  the  Fourier-Bessel  expansion  of  a  possible  model  incident  electric  field  defined  by 
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Fx(r)  =  E°H(a-r)  0  <r<a  z  =  0  (3.6) 

where  H(a  — r )  denotes  the  Heaviside  step  function.  This  excitation  source  corresponds  to  a  uniform 
field  over  a  small  circular  aperture  of  radius  r  =  a  .  The  model  could  be  used  to  represent  the  axial 
electric  field  in  a  cylindrical  resonator  operating  in  a  TM0i  mode  and  with  two  small  holes  of  radius  a  on 
the  end  caps  As  the  number  of  terms  increases,  the  Fourier-Bessel  series  converges  in  the  mean  to  the 
incident  field.  Each  Fourier  component  in  the  expansion  corresponds  to  a  particular  eigenmode  of  the 
waveguide. 


1.4 
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Figure  7  Fourier-Bessel  series  representation  of  a  uniform  incident  field  at  z=  0.  The  model  could  represent 
the  axial  electric  field  in  a  TM0i  circular  resonator  with  a  hole  of  radius  a. 


Radius  (cm) 


Substitution  of  equation  (3.4)  into  equation  (3.3)  yields 


Z 

n= 1 


( p  'y 

2  ^ 

— 

+  k2z 

K  b  ) 

y 

9  k  1 

An{t,k:)  +  ^S(cot)Yn--A:{t,kz) 


n 


P  y  OO  P  r 

<3. 

b  r?  b 


The  prime  notation  denotes  time  differentiation.  Note  that  the  conductivity  cr(V)  is  an  explicit  piecewise 
continuous  functio: 

eigen  functions  J0 


continuous  function  of  the  radial  position  r  within  the  waveguide.  If  the  conductivity  were  zero,  the 
(  p  r\ 

would  be  orthogonal  over  the  interval  [0,6]  by  Sturm-Liouville  theory. 


~  0/7 

V  b  j 
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If  we  set  the  conductivity  to  zero,  the  right-hand-side  of  equation  (3.7)  becomes  zero  and  the 
transient  model  would  reduce  to  an  infinite  set  of  linearly  independent  driven  harmonic  oscillator 
equations  (3.8)  each  corresponding  to  a  particular  waveguide  mode  and  subject  to  the  initial  conditions 
that  An(0,kz)  =  0  and  A'n(0,kz)  =  0 


2  k„c2 


<(t,  L ) + a]  ( k._ )  a,  (r,  k._ )  =  r„s(«)  o;  (k, )  = 

n 


On 

V  b  j 


+  ki 


c-  >0 


(3.8) 


J 


For  an  arbitrary  forcing  function  Si^cot ) ,  equation  (3.8)  is  readily  solved  by  taking  the  Laplace  transform 
over  time.  The  Laplace  transform  t  — >  s  is  defined  by 


A(s,kz)  =  ^A(t,k.)ex$\-st\dt  =L^A{t,kz)\  (3.9) 

o 

Taking  the  Laplace  transform  t  — >  s  of  equation  (3.8)  and  solving  for  An  ( s ,  Ay )  yields  the  solution  in  the 
[5,  A:.]  space 


A(S’K) 


2kS„  c2  S(s/co ) 
n  52+Qj(A:_)  co 


(3.10) 


The  inverse  Fourier  Sine  transform  k„  — >  z  of  equation  (3.10)  recovers  the  z-dependence. 
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\2  /  „  \2y 
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(  p 
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VU 
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=  r  „exp 


(3.11) 


Equation  (3. 1 1)  is  written  as  the  product  of  two  Laplace  transforms.  One  factor  S  |  y(0J  —  is  the  Laplace 

transform  of  the  excitation  source.  The  exponential  factor  can  be  associated  with  the  resolvent  or  Green’s 
function  for  Eqn.  (3.8).  The  Laplace  inversion  ( L  1 )  of  the  exponential  function  in  (3. 1 1 )  back  to  the 
time  domain  is  accomplished  by  the  following  substitution20 
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(3.12) 


This  Laplace  inversion  provides  a  time-dependent  Green’s  function  G{t,  z )  for  the  n-th  mode.  The  mode 
evolution  for  an  arbitrary  driving  function  can  be  obtained  by  the  Laplace  transform  convolution  theorem 


An(t,z)  =  nG{r,z)s{Kco[t -T^jdT  (3.13) 

o 

Equation  (3.13)  may  be  readily  implemented  in  MathCAD  or  Mathematica  to  explicitly  calculate  the 
mode  evolution  for  an  arbitrary  temporal  excitation  function  S  (  cot ) . 

If  we  assume  S  [t)  =  Sin  [cot) ,  equation  (3.10)  for  An  ( s,k z)  takes  the  following  form 
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(3.14) 


This  expression  can  be  inverted  back  to  the  time  domain  s  — >  t  by  the  inverse  Laplace  transform  to  give 
two  terms  for  An  ( t,k _ ) 


An(t,k:) 


co Sin  (Q  [k, )  t)  -  Cl[k_ )  Sin  [cot) 
v  kl[kz)[co2-kl2[k:)]  y 


(3.15) 


In  order  to  recover  the  z-dependence,  the  inverse  Fourier  Sine  transform  must  be  applied  to  equation 
(3.15).  The  equation  is  first  rearranged  as  follows 
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The  inversion  k_  — >  z  of  the  last  term  depends  upon  the  sign  of 
solutions: 
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The  first  equation  of  (3.17)  describes  an  evanescent  wave  that  decays  exponentially  in  the  z-direction 
from  the  source.  This  occurs  when  the  excitation  frequency  co  is  below  the  mode  cutoff  frequency.  The 
second  equation  of  (3.17)  describes  traveling  waves  excited  by  a  source  whose  frequency  is  above  the 
mode  cutoff  frequency.  The  mode  dispersion  relation  is  the  same  as  that  given  by  equation  (2.4).  The 
inversion  of  the  first  term  in  equation  (3.16)  is  under  investigation. 


We  will  now  incorporate  the  effect  of  finite  conductivity  upon  the  mode  propagation  and  loss. 
The  excitation  source  is  assumed  to  be  of  the  form  Sin  ( cot ) .  Multiplying  both  sides  of  equation  (3.7)  by 


Jn 


P0mr 

b  J 


and  noting  that  the  eigenfunctions  on  the  left-hand-side  are  orthonormal  over  the  interval 


[0,6]  yields 


~j  A.  (*,  K )  +  An  (6  K )  =  ^r1-  Sin(o>t )  -  X  im  I  aAA  |  n)  A,  (t,  K ) -jr 
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22=1 


(3.18) 
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a{r,t)\n)  =  \jS^Ai 
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f  P,.,A 


dr 


V  o  j 


Equation  (3.18)  shows  that  the  first  effect  of  the  finite  conductivity  is  to  couple  the  m-th  mode  with  all 
other  waveguide  modes.  If  the  conductivity  were  a  constant  defined  over  the  interval  [0,6] ,  the  infinite 

n-sum  on  the  right-hand-side  of  (3.18)  would  reduce  to  a  single  diagonal  term  in  which  n  =  m  and  we 
would  obtain  an  independent  damped  harmonic  oscillator  equation  analogous  to  equation  (3.8)  for  each 
mode. 
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An  approximate  solution  for  the  damping  of  the  m-th  mode  can  be  obtained  by  iteration  of 
equation  (3.18).  This  equation  can  be  rearranged  as  follows 


A'L  {t,  kz )  +  2 ymA'm  (/,  kz )  +  Q 2  [k,  )Am  (/,  kz )  =  2k?mC  ^(^0  _c- ^  J 


m  a 


Mb 


n 


m^n=\ 


N 


4(t,k:) 


(3. 


_  PvP 


Ml 


m 


The  approximate  solution  for  the  m-th  damped  mode  is  constructed  by  neglecting  the  mode  coupling 
terms  in  the  sum  on  the  right-hand-side  and  solving  the  driven  damped  oscillator  equation  by  taking  the 
Laplace  transform  over  time. 


Am(s,K)=2KTmCl 
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n  s2  +  2yms  +  Q2  (kz ) 
2  kT...c2  1 


r  co  ^ 
\s2  +  C02  J 
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is  +  /m)  +Ql{k:)~K 


k  co  ^ 
ys2+co2J 


(3.20) 


The  z-dependence  of  Am  (s,  z)  is  recovered  by  taking  the  inverse  Fourier  Sine  transform  k..  — »  z  of 
equation  (3.20) 
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The  inversion  of  equation  (3.21)  back  to  the  time  domain  is  accomplished  by  the  following  substitution 
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The  inverse  Laplace  transform  of  equation  (3.22)  may  be  written  as  a  product  of  two  time-dependent 
functions  using  the  Laplace  transform  substitution  theorem21 


k  '  [  f{s  +  /)]  =  exp(-yt)f(t) 


(3.23) 
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Using  relation  (3.23)  and  Churchill’s  Tables  of  Laplace  transforms20,  the  time-domain  solution  is  given 
by 


exp[-yj] 
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l  c) 

(OZ 


cAt  - 


UJ 


co.„At~  - 


(-] 

UJ 


(3.24) 


Equation  (3.24)  provides  the  Green’s  function  or  response  function  associated  with  the  m-th  mode.  The 
temporal  evolution  of  the  mode  coefficient  Am  ( t,z )  can  be  calculated  from  the  Laplace  convolution 

integral  for  the  sinusoidal  input  , S(t)  =  Sin  {(Ot) 
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(3.25) 


Since  the  excitation  source  S((Ot)  is  a  time-dependent  sine  function,  the  Dirac  delta  function  in  equation 
(3.25)  reduces  the  integration  over  the  first  term  to  a  damped  harmonic  oscillator  solution 


r  mexp 


r  zi 
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ry 
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(3.26) 


This  term  represents  a  spatially  damped  wave  of  initial  amplitude  Tm  traveling  at  the  speed  of  light 

within  the  plasma  tube  in  the  positive  z-direction.  This  term  corresponds  to  a  standard  perturbation 
treatment  in  which  the  attenuation  due  to  wall  conductivity  for  a  travelling  sine  wave  as  an  initial 
condition  is  calculated.  In  the  coupled  metal-plasma  tube  model,  the  amplitude  would  be  determined 
from  the  steady-state  mode  at  the  aperture  z  =  0 .  The  second  term  in  (3.25)  has  both  a  decay  and 

oscillatory  behavior  associated  with  the  Bessel  functionJj .  The  sinusoidal  driving  term  Sin  \jo(t  —  r)J 

can  be  written  as  two  terms 


Sin  rn (7  -  r ) ]  =  Sin  [ cot ] Cos  [ co z  -  Cos  [rn/  Sin[ojx 


(3.27) 
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Substitution  of  Eqn.  (3.27)  into  Eqn.  (3.25)  gives  the  following  solution 
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This  solution  is  difficult  to  interpret.  However,  we  can  consider  some  limiting  cases.  For  z=0  and  all  time 
t  >  0,  equation  (3.28)  reduces  to  rm  Sin  [cot]  which  corresponds  to  the  specified  boundary  condition  at 
z=0.  In  the  limit  oH  — »  oo  ,  the  integrals  go  to  zero  and  the  solution  becomes  a  spatially  damped  traveling 


wave  of  the  form  rmexp 


(-y„A 

Sin 

(  A 

cot  — 

V  c  ) 

l  c) 

For  intermediate  values  of  z  and  t,  the  signal  propagates 


r 


at  the  speed  of  light  c.  The  Heaviside  function  H 


T 


introduces  temporal  retardation  into  the  signal 


V  cJ 

propagation.  The  amplitude  seen  by  a  field  detector  or  probe  at  some  intermediate  position  and  time  is 
composed  of  sinusoidal  oscillations  at  frequency  CO  modulated  by  the  Bessel  function  integrals.  For 
specified  values  of  z  and  t,  these  integrals  can  be  evaluated  numerically.  Figure  8  shows  the  time 


dependence  of  the  Green’s  function  (integrand)  for  the  7C0S  (z,t)  component.  Figure  9  depicts  the  integral 


of  this  component  and  its  contribution  to  the  signal  amplitude  at  the  same  position.  The  transient 
contribution  to  the  signal  is  on  the  order  of  1 0- 1 5  %. 
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Figure  8  Green  Funetion/cos  component  at  z=50  cm.  and  f=1.0  GFIz 


Figure  9  Integral  contribution  to  signal  amplitude  at  z=50  cm.  and  1=1.0  GHz. 
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IV.  FINITE-ELEMENT  COMPUTATIONS 

In  order  to  validate  the  analytic  model  described  here,  some  preliminary  numerical  simulations 
were  performed  using  the  finite-element  code  Comsol  3.5.  The  geometry  of  the  model  is  depicted  in 
Figure  10.  The  model  geometry  assumes  a  perfectly-conducting  metal  tube  (PEC)  of  radius  5  mm  is 
connected  to  a  plasma  tube  of  the  same  interior  radius  with  a  wall  thickness  .5  mm.  Specified  TM  and  TE 
modes  above  the  tube  cut-off  frequencies  can  be  launched  from  the  excitation  port  and  propagated  into 
the  plasma  region.  The  attenuation  is  determined  by  numerically  integrating  the  axial  Poynting  vector, 
which  describes  the  power  flow,  over  the  waveguide  cross  section  at  the  plasma  inlet  port  and  at  a  second 
station  located  5  cm  from  the  end  of  the  10  cm  plasma  tube.  An  attenuation  coefficient  a  (1/m)  is 
derived  from  the  integrated  power  ratio  and  the  assumption  of  an  exponential  decay  over  the  propagation 
distance  of  5.0  cm.  The  preliminary  calculations  were  done  assuming  steady-state  propagation  in  order  to 
determine  the  finite-element  size,  grid  and  memory  requirements  for  numerical  convergence. 


PEC 


Plasma 


Mode  Excitation 
Port 


Plasma  Inlet 
Port 


Figure  10  Geometry  of  numerical  model. 


In  order  to  confine  the  wave  within  the  50  S/m  plasma  region,  the  excitation  frequencies  were  selected  by 
requiring  the  skin  depths  be  less  than  0.5  mm.  This  plasma  wall  thickness  would  correspond  roughly  to 
three  closely  packed  circular  arrays  of  100-150  pm  laser  induced  filaments.  The  results  of  these 
calculations  are  given  in  Table  III  and  Figure  11. 

Table  III  Attenuation  Calculations  TMW  Mode 


Excitation 
Frequency  (GHz) 

Skin  Depth  (mm) 
a  =  50  S/m 

Numerical  a  (1/m) 

Standard 
Perturbation  a 
(1/m) 

a  (1/m) 

Equation  (3.19) 

25 

.45 

45.6 

59.5 

53.9 

30 

.41 

42.1 

40.1 

53.6 

35 

.38 

44.2 

37.0 

53.2 

40 

.36 

47.6 

36.4 

52.6 

45 

.34 

51.2 

36.8 

51.3 

50 

.32 

55.0 

37.6 

50.6 

55 

.30 

58.5 

38.5 

49.9 

60 

.29 

62.4 

39.6 

49.2 

The  standard  perturbation  theory  underestimates  the  attenuation  in  the  plasma  waveguide.  This  is 
not  unexpected,  since  the  underlying  assumption  for  applying  the  perturbation  model  is  that  the  skin 
depth  is  small  in  comparison  with  the  plasma  wall  thickness.  The  numerical  power  integration  shows  an 
increase  in  attenuation  with  frequency.  The  loss  is  on  the  order  of  360  to  550  db/m.  Experiments 
conducted  with  this  geometry  and  frequency  range  would  not  be  able  to  detect  a  propagating  signal. 
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Perturbation  theory  indicates  that  the  loss  is  inversely  proportional  to  the  waveguide  radius.  Figure  12 
shows  the  calculated  attenuation  for  a  TEn  mode  in  a  6.0  cm  diameter  plasma  waveguide  with  a  2  mm 
thick  wall.  Under  these  conditions  the  perturbation  theory  and  the  numerical  power  integration  are  in 
reasonable  agreement  and  the  attenuation  is  reduced  by  an  order  of  magnitude. 


Figure  12  Calculated  steady-state  attenuation  in  a  plasma  (50  S/m)  waveguide  of  radius  3  cm. 
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V.  CONCLUSION 

An  analytic  model  to  describe  the  transient  propagation  of  an  electromagnetic  wave  inside  a 
hollow  constant  conductivity  plasma  annular  tube  has  been  formulated  utilizing  integral  transform 
techniques.  The  model  provides  insight  into  the  wave  attenuation,  the  nature  of  the  evanescent  and 
propagating  modes  (Eqn.  3.17),  the  dispersion  relations  and  waveguide  excitation.  The  key  result  of  the 
time-dependent  analysis  is  equation  (3.24)  which  provides  a  Green’s  function  for  the  solution  of  the 
interior  propagation  problem.  The  Green’s  function  reveals  that  the  transient  signals  seen  by  a  field 
detector  have  oscillatory  components  of  the  same  frequency  as  the  excitation  source,  but  with  time- 
dependent  modulated  amplitudes  (e.g.  Fig.  9)  that  depend  upon  the  distance  from  the  excitation  source. 
The  Green’s  function  is  a  complex  entity  and  how  it  relates  to  what  we  measure  in  any  particular 
experiment  is  under  investigation.  Although  the  Green’s  function  developed  here  is  only  applicable  to  the 
specified  geometry  of  the  cylindrical  waveguide  model,  the  integral  techniques  utilized  in  the 
development  can  be  extended  to  other  geometries.  We  propose  to  use  similar  techniques  for  the  open 
waveguide  to  understand  launching  and  possible  surface  wave  formation  on  the  laser-generated  plasma 
filaments. 

Some  preliminary  numerical  calculations  on  closed  plasma  waveguides  have  been  performed 
utilizing  a  finite-element  code.  The  simplest  and  probably  the  best  possible  configuration  for  guided 
wave  propagation  is  a  steady-state.  Published  papers2'3  have  utilized  conventional  perturbation  theory, 
which  assumes  a  given  mode  is  propagating  within  the  waveguide,  to  predict  the  attenuation  of  that  mode 
within  the  hollow  plasma  tube.  We  have  compared  both  the  perturbation  theory  and  the  finite-element 
calculations  for  wave  attenuation  in  a  50  S/m  plasma.  In  order  to  confine  the  wave  within  the  waveguide 
over  the  4-10  GHz  frequency,  the  annular  plasma  thickness  should  be  on  the  order  of  a  few  millimeters 
with  a  waveguide  radius  of  3.0  cm.  The  calculated  attenuation  in  this  configuration  is  on  the  order  of  12- 
20  db/m  for  the  TEn  mode.  For  this  particular  geometry  and  mode,  the  perturbation  calculations  and  the 
numerical  integrations  agree  reasonably  well  as  shown  in  Figure  12.  As  the  waveguide  diameter  is 
reduced,  the  attenuation  increases.  At  the  higher  frequencies  and  for  smaller  diameter  waveguide,  the 
perturbation  theory  and  the  numerical  computations  for  the  attenuation  are  in  poor  agreement  as  shown  in 
Figure  11. 

Finally,  we  note  that  an  alternative  approach  to  guiding  the  electromagnetic  wave  is  to  use  a 
waveguide  whose  radius  is  much  larger  than  the  wavelength  of  interest.  The  guiding  mechanism  is  based 
on  the  effect  of  total  reflection  at  the  plasma-air  interface.  Signals  at  38  GHz  have  been  directed  over 
distance  of  60  m  by  this  approach  which  is  referred  to  as  a  sliding  mode2’3. 

Part  of  the  motivation  for  doing  the  steady-state  calculations  and  comparing  them  with 
perturbation  theory  is  to  confirm  the  accuracy  and  predictive  character  of  the  numerical  model.  The 
annular  tube  waveguide  problem  has  both  closed-form  and  perturbation  solutions  which  can  be  compared 
with  the  numerical  results  to  determine  the  appropriate  input  conditions  such  as  grid  size,  plasma  wall 
thickness,  propagation  length  and  outer  region  boundary  conditions.  If  the  annular  plasma  wall  were 
replaced  with  a  ring  of  plasma  filaments,  there  are  no  analytic  expressions  for  the  solutions  to  compare 
the  numerical  predictions  with.  The  problem  is  also  complicated  by  the  fact  that  the  computational 
domain  is  now  open  to  free  space  and  the  wave  can  “leak”  out  into  the  surroundings.  Thus,  one  of  our 
first  objectives  was  to  establish  the  conditions  under  which  the  numerical  work  agrees  with  a  known 
analytic  model  and  where  it  becomes  invalid. 
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